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Abstract. Complete next-to-leading (chiral) order (NLO) expressions for the strong 
isospin-breaking (IB) contributions in K ^ tttt are used to discuss (1) for CP-even, 
the impact on the magnitude of the A/ =1/2 Rule, and (2) for CP-odd, the strong IB 
correction, flgt, for the gluonic penguin contribution to e'/e, with particular emphasis 
on the strong low-energy constant (LEC) and loop contributions, numerical values for 
which are model-independent at NLO. 



In the presence of IB, the standard isospin decomposition of the tt+tt 
TT+TT^, 7r°7r° decay amphtudes, A^q, A+_ and Aqq, becomes [1] 







Aoo = Vl/3Aoe^*° - [v/2/3]A2e**^ 
= /l/SAoe^*" + [l/v^]A2e**^ 



A+o = [VsmA'.e"''^ (1) 

Ignoring AI = 2 electromagnetic (EM) channel coupling, $/ are the tttt phases. In 
general, EM- and strong-IB- induced A J = 5/2 effects make [Agl 7^ IA2I. Aq, A2 
can be chosen real in the absence of CP violation. 

Since \Aq\ ~ 22|A2|, IB "leakage" of the large octet (/ = 0) amplitude into the 
AI = 3/2 (/ = 2) amplitude can be numerically significant. EM leakage has been 
computed to NLO in Ref. [1] (see also [2]); we compute the NLO strong leakage for 
both CP-even and gluonic-penguin-mediated CP-odd transitions. Strong cancel- 
lation between gluonic penguin [Oq) and electroweak penguin contributions make 
Standard Model predictions for e'/e sensitive to the latter [3]. 

At leading (chiral) order (LO), the ratio, 6A2/AQ, with 6A2 the strong octet 
leakage contribution to A2, is unambiguous. The IC part of A2, ^2*^, is decreased 
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TABLE 1. Strong octet and EM IB leakage contributions in units of 10"^ MeV. 
The isospin-conserving (IC) and LO IB fits yield A2 = A2 = -2.1 x 10~^ MeV and 
-2.4 X 10~^ MeV, respectively. 





Source 










(8) 

(EM) 


(-1.56 ±0.63) + 
(-1.27 ±0.40) - 


(0.42±0.05)i 
(1.28±0.02)i 


(-1.56 ±0.63) - 
(0.70 ±0.73) - 


f (0.42±0.05)i 
(0.07±0.04)i 



by = 13% once strong IB is included. The LO Oq suppression in e'/e is 

1 — Qst- Recent analyses of e'/e employ ilgt — 0.25 ±0.08 [3]. The difference reflects 
estimates of NLO effects mediated by an intermediate rj' through the K'^ — > tt'^t)' 
transition [4]. This effect, in conventional ChPT (involving only the vr, K and rj as 
explicit degrees of freedom), represents one contribution to the CP-even (CP-odd) 
NLO weak LEC's, {Ej^) [5], but does not exhaust such contributions. 

Let us flx notation. At LO, the low-energy representations of the dominant 
CP-even octet operator and CP-odd (— ) gluonic penguin operator are [5] 

(2) 

where = \q, = A7, and U = exp {iX • tt), as usual. A complete set of NLO 
contributions is obtained by evaluating the sum of (i) graphs with a single LO 
weak vertex (proportional to c^) and one of the external legs dressed by a single 
strong NLO vertex (proportional to one of the NLO strong LEC's, Lk) ("strong 
LEC contributions"); (ii) 1-loop graphs involving a single LO weak vertex and 
strong vertices only of LO ("loop contributions"); and (iii) tree graphs involving 
a single NLO weak vertex [5] ("NLO weak LEC contributions"). Only the sum of 
all three classes of contribution is renormalization-scale-independent and physically 
meaningful. Prom Eq. (2), it is immediately obvious that LO CP-even and CP-odd 
decay vertices are related through the substitution c+ ^ ic~. The loop-plus- 
strong-LEC part of the NLO contribution to the ratio /Aq is, therefore, the same 
for the CP-even and CP-odd cases. Because, moreover, the LO LEC's, c^, cancel 
in the ratio, this contribution is, though scale-dependent, modeZ-independent. 

The strong LEC contributions to 5A2/AQ have recently been discussed in Ref. [6]. 
They include an ?7'-mcdiated n^-r] mixing term proportional to Lj which is, however, 
empirically, almost completely cancelled by an accompanying Lg contribution [6]. 
The loop contributions have been evaluated in Ref. [7]. Since the full set of {E^) 
cannot be determined empirically at present, one is forced to use model values. The 
models used arc discussed below. 

For the CP-cvcn case, the leakage contributions, 6A2 = SA2, are given in Table 
1 [7]. The corresponding EM contributions [1] are shown for comparison. The 
errors reflect uncertainties in the estimates of the NLO LEC's (see Ref. [7] for 
details, and references to the models employed for the E'^). Denoting the ratio 
of fitted LO 27-plet to octet weak LEC's obtained neglecting, or including, IB by 
ric, or riB, respectively, we find Rib = tib/tic = 0.963 ± 0.029 ± 0.010 ± 0.034. 



The errors correspond, respectively, to model dependence of the E'j^, uncertainties 
in -Bo(m(i — to„), and the EM uncertainties of Refs. [1]. The deviation from 1 
is significantly smaller than at LO (where Rjb = 0.870). The true AJ = 1/2 
rule enhancement, tjb, can thus be approximated by rjc to an accuracy of better 
than 10%. The A/ = 5/2 contribution (dominantly EM in character [7]), leads to 
^2/^2 — 1-094 ± 0.039 7^ 1, and significantly exacerbates the phase discrepancy 
problem for the neutral K decays [7]. 

For the CP-odd case, Clgt = ujlra5A2/ ImA^ (u = ReAo/ReA2 ~ 22). ImAo is 
associated with Oq. At LO, [flst]LO = At NLO 



St 



stlLO 



= [1 + i?2 - Ro] ■ (3) 



The superscript {NLO; ND) indicates the sum of non-dispersion NLO contribu- 
tions (involving NLO weak and strong LEC's and the non-dispersive parts of loop 
graphs). Neither the NLO J = isospin-conserving nor NLO 7 = 2 IB leakage 
combinations are known. The NLO dispersive contributions create phases con- 
sistent with Watson's theorem. Although the positive 7 = phases correspond 
to attractive FSI, the Ej^ terms may, nonetheless, make Im^o smaller at NLO 
than the LO (see comments on Ref. [8] in Ref. [9] for a related discussion). If, 
however, NLO effects do enhance Im Aq (decreasing the level of Oq-Os cancellation 
and increasing e'/e) Qgt will be simultaneously suppressed, further amplifying this 
increase. The known NLO contributions (loops and strong LEG terms) give con- 
tributions -0.24(-0.31) to i?2 and -0.02(+0.42) to Rq, at scale n = m^(mp) [10]. 
For the (model-dependent) E'j^ contributions, we considered the possibility of us- 
ing 3 models: the weak deformation model (WDM) [11], the chiral quark model 
(ChQM) [12], and the scalar saturation model of Ref. [13] (SSM). Ref. [13] was 
the first to argue explicitly that NLO weak LEG contributions other than those 
mediated by the r)' might be numerically important. A large negative contribution 
to i?2 — -Ro was obtained which, if correct, would significantly enhance the Oq con- 
tribution to e'/e. In contrast, the two earlier models produce positive, or near-zero, 
E^ contributions to R2 — Rq- Since there is a problem with the SSM estimate^ we 
have used only the WDM and GhQM. Taking the WDM model E^ to correspond 



^) The model values, of the are obtained by (1) approximating the renormalized 4-quark 
operator 0% as the product of renormalized scalar densities (the factorization approximation) 
and (2) dropping (divergent) seagull terms in the low-energy representation of this product. 

Two types of contribution to result, E^^ (proportional to the strong 6*'' order LEC's), 

and E^ (proportional to a product of 2 4*'' order strong LEC's). As an example, E^ = 

IdV^G FV*gVud Im {C())Bq [13]. In Ref. [13] it is assumed that the finite part of this expression 
is proportional to the square, [ig]^' the finite part of Ls- This, however, is incorrect. In 
dimensional regularization, the Laurent expansion of Lg begins at 0[l/{d— 4)]; the finite part of 
Ll is thus actually [L^]'^ + 5Lg~^^ /SSAtt'^ , with Lg~^^ the coefficient of the 0[rf-4] term in the Lg 
Laurent expansion. ig~^^ enters physical processes beginning at 6*^ order in the chiral expansion, 



to a scale n ~ rup, we find, combining with the strong LEG and loop contributions 
given above, 1 + i?2 - -Ro = 0.27 in the WDM and 0.64 in the ChQM. Taking the 

more microscopic ChQM to provide a central value, and the difference between 
the WDM and ChQM to provide a minimal measure of model dependence in the 
theoretical result, we obtain, at NLO, Qgt = 0.08 ± 0.05, significantly smaller than 
both the conventionally employed value 0.25 ±0.08, and that, 0.16 ±0.03, obtained 
in the estimate of Ref. [6], based on the strong LEC contributions only. The central 
value above, combined with conventional central values for the 5-factors, leads to 
a ~ 50% increase in the predicted value for e'/e. To be conservative, we would 
propose using this lower central value with an even larger error estimate, in all 
future calculations of Standard Model values for e'/ e. 
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and hence is on the same footing as the other 6*'' order strong LEC's appearing in E,^ and 

retained in Ref. [13]. Since, it is inconsistent to drop the L'j^^\ and the SSM provides no means 
of estimating their values, the numerical estimates given in Ref. [13] are incomplete, and cannot 
be used. 



